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Abstract. - We investigate the role of the form of the spatial diffusion coefficient in shock ac- 
celeration of fast particles. Referring to non-classical diffusion and using the results of numerical 
(hybrid) simulations tailored for the downstream shock population in quasi-perpendicualr high- 
Mach number collisionless shocks to which we apply the theory, we demonstrate that the inferred 
diffusion coefficients are in excellent agreement with the requirements of the theory and its pre- 
dictions. Diffusion in the collisionless regime turns out to be non-classical Gibbsian (Levy flight), 
time-dependent though weak. 
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J Introduction. — Diffusive shock acceleration theory of 
particles out of the thermal flow to high energies makes ex- 
tensive use of spatial diffusion coefficients, assuming that 
these coefficients are generated self-consistently in the in- 
teraction of plasma waves and particles. The multiple tran- 
sitions of a fast particle back and forth across the shock 
under the assumption that the scattering is a stochastic 
process suggests a diffusive approach to the acceleration 
mechanism. In the simple model where the waves prop- 
agate just parallel and antiparallel to the magnetic field 
this diffusion tensor is given by 
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where 



der the assumption that v^/u), 
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(1) 



j. These expressions hold un- 
<C 1, withj/^ = 2D^^/(1- 



fx") the pitch-angle scattering frequency, fi = cos a, I?^^ 
the pitch angle diffusion coefficient, and a the pitch an- 
gle. Then <C 1 as well. Usually the perpendicular 
diffusion coefficient k± is also neglected. However, while 
the former can be justified somehow, k± is implicitly de- 



Present address: International Space Science Institute, Bern, 
Switzerland, e-mail: treumaiin@issibern.ch 



termined [3] by an expression like 

■■2 5(k)d3k 



Kj_ 
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K± + k?,K\\ + v/Xb\\ + 7(k) 



(2) 



where a is a constant of proportionahty that is determined 
from simulations, A(,|| is the parallel correlation length 
of the fluctuating magnetic fields, and S{k) is defined 
through S{li,t') = 5'(k)r(k,f'), with T = exp[-7(k)i'], 
being the spatial correlation of the transverse magnetic 
fluctuations 

(6,[x(0),0]6,[x(O,i']) = j R{y,t')P{y\t')dy (3) 

R{y,t') is the two-point, two-time correlation function, 

and P{y\t') is the probability density for a particle at time 
t' being displaced by the amount y in transverse direction. 
Moreover, the two-time parallel-velocity autocorrelation 
is modelled by the isotropic assumption {vi\{0)v^^{t')) = 
{v"^ /3) exp{—vt' / Xb\\) , and it is assumed that 7 = 0. With 
P being a symmetric Gaussian distribution for the par- 
ticle trajectory (assuming that the displacement is at all 
times diffusive) , the average over the exponential wave fac- 
tor (exp[ik • x{t')]) which is introduced through the mag- 
netic field fluctuations, yields just the cumulative average 
(exp[ik • x(i')]) = exjp[{—k'j_K± — A:jj«;||)t']. 

This Non-Linear Gyro-Center (NLGC) diffusion model 
is close to reality as it takes account not only of the random 
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walk of the fluctuating magnetic field lines but also of the 
distortion of the particle orbits in this random walk. Dif- 
fusion of particles is caused by decorrelation of their orbits 
from the fluctuating magnetic field. A further development 
in this theory [9] does not use just Alfvenic turbulence 
but takes into account also the contribution of compres- 
sive modes, but the deviations from the NLGC model and 
the improvement over it are small and can mostly be ne- 
glected. 

The agreement between the NLGC diffusion model and 
the two-dimensional numerical simulation based on two- 
dimensional magnetic field turbulence is reasonable when 
compared with other models. However, the perpendicular 
diffusion coefficient is considerably less than quasilinear 
theory predicts, and the deviations from quasilinear the- 
ory become susceptible already at rd ^ 0.1 A. Still this 
model is not yet self-consisten as it ignores the feedback 
of the particles on the magnetic field turbulence which - 
at least to some degree - is a function of the presence of 
the energetic particles, as suggested by self-consistent the- 
ory [1]. Determination of the diffusion coefficient at this 
stage is a necessary intermediate step. A complete solu- 
tion of the problem including parallel and perpendicular 
diffusion for the relevant waves as well as energy diffu- 
sion can be expected to come only from three-dimensional 
full-particle PIC simulations of shock formation including 
particle diffusion and acceleration to high energies. 

Non-classical difFusion ('super-diffusion'). The 

above diffusion coefficients are all based on the assump- 
tion of classical diffusion. Since shocks are narrow tran- 
sitions from one plasma state to another one particle in- 
teractions in the vicinity of shocks might be subject to 
statistics of extremes and not to classical statistics. In 
this case the diffusion becomes time-dependent as no fi- 
nal state is reached in the process of particle scattering. 
The probability P{y\t') — > P{y\t'; u) is not anymore Gaus- 
sian but has a long tail of power —i/ that extends toward 
the rare large excursions and long waiting times between 
the excursions. The theory of such distributions goes back 
to [2]. More contemporary developments based on frac- 
tal theory have been given in [6] and [4], and a statistical 
mechanical argument has been developed in [8]. An early 
application of generalised Gibbsian statistical mechanics 
(respectively Levy flight theory) to the derivation of the 
parallel and perpendicular diffusion coefficients [7] yields 
a time dependent (parallel) mean square particle displace- 
ment {z(0)z{t)) = K|| {d, t; v)t which leads to a time depen- 
dent parallel (or non-magnetised) particle diffusion coeffi- 
cient 

K\\{d,t;v) = K||aK(l/-rf/2)](l/ant)" (4) 

a = (4z/-2rf-l)-^ (5) 

Here d is the dimension of the system, i^an is an anomalous 
collision frequency that is usually nonzero, and Kyci is the 
classical diffusivity based on i/an- The condition that P is a 
positive valued probability distribution is that v—djl > 1. 



The diffusion coefficient then scales as /ty oc Tii'^n^t", 
which for the scaling with the collision frequency gives a 
scaling exponent | < |a— 1| < 1- In order to find the 
scaling of the perpendicular diffusion coefficient one can 
use the ordinary classical formula 

K_L = K||/(l + W^Jl/an^) (6) 

(a more precise theory should usc^ tlw^ integral equation for 
K± given in Eq. (2) instead of this classical approximation) 
which yields 

= Jj;. (J^ant) (7) 

U - d/2 \uJci ) TUiUg^n 

Under the condition that fc'an oc ojd this yields the scal- 
ing K± OC KBioJcit)" with < a < I and kb the Bohm 
diffusion coefficient. 

It is very important to realise that the entire physics of 
deviation from purely stochastic diffusion processes is con- 
tained in this surprisingly extraordinarily narrow range of 
exponents. Very precise determination of this exponent is 
therefore crucial for elucidating the actual physics of dif- 
fusion and hence also diffusive particle acceleration. Be- 
low we will discover in numerical simulations of energetic 
particle diffusion near shocks that the non-classical limit 
might indeed be realised there. 

This kind of transverse diffusion increases with time 
and, thus, the displacement of the particles grows at a 
faster rate than any classical transverse displacement and 
even faster than under Bohm diffusion kb- However, if the 
system would have sufficient time for reaching a final state 
it would end up in the limit of classical diffusion with lin- 
ear growing root mean square displacements. This means 
that under non-classical conditions the system evolves dif- 
fusively faster than it would evolve under classical condi- 
tions, while given sufficient time it would settle in classical 
diffusion. A diffusion of this kind is called 'super-diffusion'. 
In an infinitely extended system it will be realised only 
temporarily in the initial state for times shorter than the 
typical classical collision time. We will return to this im- 
portant statement in the discussion section. 

In a spatially limited system like the shock transition 
the final stationary diffusion state might not be reached, 
however, since the particles have not sufficient time to un- 
dergo classical collisions. In this case the above classical 
theories require correction for time dependent diffusion, 
and a stationary state is achieved only when the time 
dependent diffusion is balanced by losses at the moment 
when the mean displacement exceeds the size of the 
downstream region. From the definition of k± = A^/ro 
one may estimate the limiting diffusion time td as 




which yields the following limitations on the anomalous 
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Fig. 1: Two dimensional numerical simulation result of the mean downstream perpendicular displacement of quasi-perpendicular {Obu = 
87°) supercritical (Ma = 4) shock-accelerated ions as function of simulation time (simulation data taken from [5], courtesy American 
Geophysical Union). Left: The displacement performs an initial damped oscillation before settling in a continuous diffusive increase at 

about t 40ll)~^^ . The further evolution deviates slightly from classical (linear) increase, following a (Ax^) oc (tuid)^'^'^ power law. This 
is very close to a power of g, suggesting that the particle diffusion process is in fact 'super-diffusive' [7]. Right: Late time trajectory of 
an arbitrary ion projected into the plane perpendicular to the mean magnetic field. The ion orbit in the superposition of the ambient and 
wave magnetic field is not a smooth stochastic trajectory, ft consists of waiting (trapped gyrating) parts and parts when the ion suddenly 
jumps ahead a long distance as is typical for rare extreme event statistics. 



collision frequency: 

The left limit holds for the extreme case v — d/2 = 1, the 

right for v ^ oo (the latter corresponding to a Gaussian 
probability distribution). For the Earth's bow shock with 
the downstream region being the magnetosheath of width 
Ad ~ 2 Re the relative diffusion times will be roughly 
i'anTD ^ 2 X 10'' (here IRe = 6436 km is the Earth ra- 
dius). Measuring the limiting diffusion time td of the en- 
ergetic particles provides an opportunity to determine the 
anomalous collision frequency u^n that governs the diffu- 
sive interaction. 

Shock simulation of particle dlfTuslon. — In order to 

confirm the diffusive nature of the ion acceleration pro- 
cess in supercritical shocks, three-dimensional numerical 
hybrid simulations (with particle ions and a neutralis- 
ing electron fluid) have been performed [5] measuring the 
root-mean-squaxe high energy ion displacement (Ax^) in 
the direction perpendicular to the magnetic field in the 
downstream region as a function of simulation time. The 
simulations corresponded to the region downstream of a 
quasi-perpendicular supercritical shock of Mach number 
A4a = 4: and shock- normal angle Obu = 87°. Due to re- 
strictions of computing power the diffusion process could 
be followed only just up to simulations times roughly 



The result of these simulations is reproduced here in 
Figure 1 with simulation data taken from [5]. Here (Ax^) 
has been plotted as a function of time ttOci- Initially the 
displacement performs a large amplitude damped oscilla- 
tion until the diffusive equilibrium is attained at about 
tiVci ^ 40. For later times the particle displacement in- 
creases continuously. However, the exponent of the in- 
crease is found not to be unity as was expected for classical 
diffusion. It is rather larger, being 1.17, which is very close 
to g identifying the diffusion process as 'super-diffusion'. 
(Note that, because of the large number ~ 6.3 x 10^ of 
macro-particles used in the simulation of which 525000 
have high (Hi(^rgi(^s and contribute to the determination 
of the mean displacement, and because of the high time 
resolution, the statistical error of the measurement is less 
than the line width!) In [5] this discrepancy is well noted 
while, nevertheless, it is insisted on interpreting the sim- 
ulations in terms of the classical diffusion picture. Here 
we use these data in order to check whether they fit the 
non-classical diffusion theory. 

Above it has been shown that the exponent of the mean 
displacement in anomalous or 'super-diffusion' is 1+a with 
a = (4i/-2d-l)~^, as given in Eq. (4) and < a < ^. The 
lower bound corresponds to classical diffusion (the linear 
increase of the displacement in the log-log representation 
of Figure 1). In the simulated case a ~ ^ lies clearly in the 
permitted range of exponents. We may use this value to- 
gether with the known dimensionality d = 3 of the simula- 
tion in determining the value of the anomalous parameter 
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Fig. 2: Evolution of the anomalous perpendicular diffusion coefficient under conditions of Levy statistics. Left: The diffusion coefficient is 
a power law increasing with time t with anomalous collision frequency i/an. Its maximum slope is +|. Approaching collision time u^j^ it 
merges into the constant classical diffusion coefficient (shaded region), it Right: Evolution of the rms perpendicular displacement (x^) as 
function of time. The maximum slope of the power law is ~ 1.3. When classical diffusion takes over the displacement evolves linearly with 
time. 



1/ = ^ = 3.25. (If using the exact value a = 0.17""'^, we 
obtain i/ = 3.22.) This value is sufficiently far above the 
marginal value of z/ = | for the three-dimensional case (or 
f = I for the two-dimensional case). For the basic theory 
see [8]. It is, however, also far enough below the classi- 
cal diffusive limit u ^ oo thus identifying the diffusion 
process as indeed being anomalous, non-stochastic (non- 
Markovian) and super-diffusive. Since these simulations 
are completely collision-free, particle diffusion is entirely 
determined by anomalous processes that are mediated by 
the self-consistently excited wave spectrum downstream of 
the shock. 

In view of this reasoning it is instructive to look at the 
trajectory an arbitrarily chosen fast ion performs in the 
magnetic fluctuation field. The right part of Figure 1 taken 
from [5] shows such an orbit in the plane perpendicular to 
the magnetic field at a late time 100 < tuid < 160, when no 
average diffusive displacement was determined anymore. 
Clearly, the trajectory is far from being a smooth diffusive 
particle orbit. It rather consists of a sequence of gyrating 
sections and sections when the particle breaks out of gy- 
ration. When performing the former sections, the particle 
is 'waiting' at its guiding centre location, approximately 
conserving its magnetic moment and performing a slightly 
modulated gyration around the magnetic field which is due 
to weak interaction with the wave field. In this regime it is 
probably interacting only very weakly with the wave fields. 
During the break-out sections it suddenly jumps out in or- 
der to occupy another waiting position. The break-out is 
most probably caused by a brief intense wave-particle in- 
teraction. (This could have but has not been check in the 
simulations of [5].) Clearly, during this break-out the mag- 
netic moment of the particle is not anymore conserved. 
The instantaneous location of the chosen particle is pro- 
gressing only in a;-direction. Other particles also break out 
into direction y. On a short time scale such a process can- 



not really be described anymore as being stochastic. 

Discussion. — The non-stochasticity is reflected in the 

time dependence of the perpendicular diffusion coefficient. 
However, as we have noted above, this time dependence 
will last only as long as the time remains to be shorter 
than the collision time, when the diffusion coefficient will 
assume its constant classical value. This reasoning implies 
that the actual diffusion in the collisionless regime is very 
slow. This has been noted (with surprise, because refer- 
ring to classical diffusion) from the simulations [5]. The 
diffusion is in fact much weaker than classical diffusion 
even though the increase of the displacement is substan- 
tially faster than under classical diffusive conditions. Using 
classical diffusion in diffusive acceleration grossly overesti- 
mates the diffusive effect. The term 'super-diffusion' does 
not apply to the strength of diffusion but just to its time- 
dependence. 

The excellent agreement between the numerical simu- 
lations and the theoretical predictions and limitations on 
the transverse diffusion coefficient in super-diffusion raises 
the question whether and for how long diffusion will re- 
main to be anomalous. Obviously the system which we 
have considered is completely collisionless. Also, the sim- 
ulations which we have referred to and used in order to 
estimate the time dependence of the anomalous diffusion 
coefficient, have been performed under strictly collisionless 
conditions. (In fact, because of computational limitations 
they could not be extended into longer times.) 

It is reasonable to assume that any real collisionless sys- 
tem will ultimately see collisions between the particles if 
only permitted to wait sufficiently long. This will hap- 
pen when t ~ Tc = where Tc is the binary colli- 
sion time, and Uc the binary collision frequency. At this 
time classical diffusion sets on, and k Kd = const, 
(a;^) = {x^{t = Tc))t. The constant of proportional- 
ity herein is the average displacement the particles have 
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reached at time t = Tc- How this displacement is reached 
is, however, not described by classical theory, as in classical 
diffusion theory the time is measured in units of the col- 
lision time Tc- Rather non-classical diffusive processes like 
the one proposed in [7] and used in this Letter are respon- 
sible for the particles to provide the initial condition for 
classical diffusion. This is graphically shown in Figure 2, 
where the evolution of the mean-square displacement and 
the corresponding evolution of the diffusion coefficient are 
shown for times t < Tc- Even though the evolution of the 
mean-square displacement with time is faster than under 
classical conditions, diffusion is weak as it reaches the ini- 
tial condition for classical diffusion only at collision time 

Tc- 

Classical diffusion is not starting from zero as usually 
believed, but rather from the value of displacement pro- 
vided by the non-classical diffusion process. Correspond- 
ingly, the il diffusion coefficient is smaller than 
the classical diffusion coefficient even though it evolves 
in time for catching up with the classical diffusion coef- 
ficient at time t = t,.. Hence, even though non-c;lassic;al 
diffusion evolves with time, the absolute value of the non- 
classical diffusion coefficient is smaller than the value of 
the classical diffusion coefficient in the same medium after 
coUisional conditions are reached. Super-diffusion is a fast 
process, but it is weaker than classical diffusion. 

We finally comment on the narrow range of powers a 
inferred both from super-diffusive non-classical diffusion 



theory and simulations. Even though the range of the re- 
alised powers is so small, it contains all of the important 
physics. In general power laws seem to be rather insensi- 
tive to the subtleties of the underlying physical processes. 
In observations this fact is generally taken not as serious. 
However, slight deviations in power may imply completely 
different physics as our example of shock accelerated dif- 
fusion has shown. Determination of powers must be done 
with care, and caution must be applied in their physical 
interpretation. 
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